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I. INTRODUCTION 


1 . 1 Prel imi nary Remarks 

This report presents a computer-based method of analysis of spur 
gear system dynamics. The report is based upon research conducted by 
the authors at the University of Cincinnati with the support of the NASA 
Lewis Research Center from 1982 to 1985 under Grant NSG 3188. The 
report represents a portion of the first author's doctoral dissertation. 

Knowledge of the dynamic effects in gear systems has been of 
increasing interest--stimulated by demands for stronger, higher-speed, 
improved-performance, and longer-lived systems. There have been 
numerous research efforts directed toward gear dynamic analysis. Still, 
the basic behavior of the gear system has not been satisfactorily 
understood. 

For example, in industrial settings, a high performance gear system 
is often obtained by overdesigning at the sacrifice of cost, material, 
and compactness. In aerospace or military applications, where weight is 
a premium, geared systems are often designed under conditions very close 
to the failure limits, thus introducing uncertainties in performance and 
life prediction. Moreover, gear systems are generally designed using 
static analyses. However, when gear systems are operating at high- 
speed, there are several factors affecting their performance which are 
enhanced by the speed. Specifically, the gear behavior can be affected 
by such items as: 

1) Torsional stiffness of the gear shafts. 

2) Gear tooth loading and deflection. 

3) Gear tooth spacing and profile errors 

4) Speed of the rotating bodies. 
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5) Mounting alignment. 

6) Dynamic balance of rotating elements. 

7) The mass of the gears and shafts. 

8. The mass of the driving (power) and driven (load) 
elements. 

There is not full agreement by researchers on the best methods for 
evaluating dynamic load effects. Hence, a gear designer is often 
confronted with conflicting theories. Therefore designers have had to 
rely on past experience, service safety factors, and upon experimental 
data with limited range of applicability. 

If major advances in mechanical reliability, optimum performance 
and design automation are to be obtained, more in-depth understanding of 
the dynamic behavior is needed. Computer analysis and experimental 
investigation are two suitable approaches for this study. Experiments 
which produce reliable results are costly involving time consuming 
procedures. Experiments are generally more expensive, both in time and 
money, than computer analysis. Hence, if reliable computer software can 
be developed, taking into account the parameters used in experiments, 
many experimental hours in the analysis and design of spur gear systems 
could be saved. 

1 .2 Literature Survey 

Research efforts on gear system dynamics have been conducted for 
many years. In 1892, Lewis [1]* recognized that the instantaneous load 


*Numbers in brackets refer to References at the end of the Report. 
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of the tooth was affected by the velocity of the system. In 192‘>, a 
large experimental program was started by a research committee, headed 
by Earl Buckingham [3], and endorsed by the American Society of 
Mechanical Engineers. They published the first authoritative report on 
gear dynamics in 1931. This report presented a procedure on determining 
the so-called dynamic load increment due to mesh dynamics and gear tooth 
errors . 

In 1959, Attia [4] performed experiment to determine actual 
instantaneous loading. He found that Buckingham's results gave more 
conservative values. 

In 1958, Niemann and Retti [5] found that larger masses caused 
higher dynamic loads, but as the average load became larger the effect 
of larger masses became less important. They also found that very 
heavily loaded gear systems showed no appreciable dynamic load 
increment, whereas in lightly and moderately loaded gear systems there 
were considerable dynamic load increments. In 1958, Harris [6] 
suggested that for gear systems isolated from external stimuli, there 
are three internal sources of dynamic loads: 1) Error in the velocity 

ratio measured under the working load; 2) Parametric excitation due to 
stiffness variation of the gear teeth; and 3) Non-1 i neari ty of tooth 
stiffness when contact is lost. In 1970, Houser and Seireg [7] developed 
a generalized dynamic factor formula for spur and helical gears 
operating away from system resonances. The formula took into 
consideration the gear geometry and manufacturing parameters as well as 
the dynamic characteristics of the system. 

In 1972, Ichimaru and Hirano [8] analyzed heavy-loaded spur gear 

systems with manufacturing errors under different operating conditions. 

They found that the change in tooth profile showed a characteristic 

trend to decrease dynamic load. In 1978, Cornell arid Westervelt [9] 
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presented a closed form solution for a dynamic model of spur gear system 
and showed that tooth profile modification, system inertia arid damping, 
and system critical speeds, can have significant effect upon the dynamic 
loads. In 1981, Kasuba and Evans [10] presented a large scale digitized 
extended gear modeling procedure to analyze spur gear systems for both 
static and dynamic conditions. Their results indicated that gear mesh 
stiffness is probably the key element in the analysis of gear train 
dynamics. They showed that the gears and the adjacent drive and load 
systems can be designed for optimum performance in terms of minimum 
allowable dynamic loads, for a wide range of operating speeds. 

In 1981, Wang and Cheng [11] developed another dynamic load 
response algorithm. They reported that the dynamic load is highly 
dependent on the operating speed. Nagaya and Uematsu [12] stated that 
because the contact point moves along the involute profile, the dynamic 
response should be considered as a function of both the position and 
speed of the moving load. In 1982, Terauchi, et al. [13] studied the 
effect of tooth profile modifications on the dynamic load of spur gear 
systems. According to their results, the dynamic load clearly decreased 
with proper profile modifications. 

1 . 3 Objective 

The objective of the research reported herein is the development of 
a theoretical basis and an associated comptuer-aided design procedure, 
for studying dynamic behavior of spur gear systems. Procedures for the 
research include the following steps: 
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1 . 


The development of a mathematical model of spur gear systems 
consisting of the gears, the shafts, and the connected masses. 

2. The evaluation of the properties of the system components. 

3. The development of governing equations for the entire system. 

4. The development of a static and dynamic analyses of the system. 

5. An examination of the effects upon system behavior of the following 

parameters : 

Mass moments of inertia 
Stiffness of components 
Operating speeds of system 
Damping factors 
System Natural Frequency 
Contact ratios 

6. The development of system component modification for design 

optimization. 

The computer software developed herein is intended to be used by 
engineer and designers in the development of actual gear systems. The 
developed software is believed to be among the most detailed and 
advanced prepared to date. 
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II. THE SYSTEM AND ITS PROPERTIES 


1 1 . 1 System Definition and Methods of Analysis 

Figure II. 1 shows a sketch and a schematic model of NASA 3000-hp 
transmi ssion stand. The system is very el aborate. However, one can 
view the system as containing three basic elements: The first is the 
gear; the second is the shaft; and the third is the mass connected to 
the shaft. Hence, an elaborate system can be modelled by components 
consisting of simple gear set. Figure II . 2 illustrates a model of a 
simple spur gear system with those three basic elements. 

The system properties such as stiffness, inertias, dampings, 
frictions, contact ratio, need to be determined to develop the dynamic 
analysis . 

The static properties of the system and its components may be 
obtained from the literature in gearing, and from the principles of 
strength of materials, mathematics, lubrication, vibration analysis, and 
finite element methods. These properties are useful in the sequel in 
conducting the dynamic analysis. 

Given the model, the governing equations of the system are written 
and integrated. It is assumed that the dynamic activity of the gears is 
confined to the rotating plane of the gears. The rotating axes are 
assumed to be symmetrical. Out-of-plane twisting and misalignment 
effects are neglected. 

A parametric study is performed to examine the relative effects of 
rotating speed, applied load, lubrication and damping, mass, stiffness, 
diametral pitch and contact ratio. 

Finally, analaytical procedures are developed for computer-aided 
design of the system. Specifically, the system is modified to obtain 
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Figure II. 2. A Simple Spur Gear System Model 
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improved performance based upon the results of the simulation. These 
are discussed in a subsequent report. 


11.2 Spur Gear Geometry: The Involute Curve and T t.s Use as a Tooth Profi le 

The Involute curve is used almost exclusively for spur gear tooth 
profiles. The involute provides numerous kinematic advantages. The 
involue generally described as the curve generated by the locus of the 
end of a line unwound from the circumference of a circle. See Figure 
II. 3. The circle from which the string is unwound is called the "base 
circle." 

The equations of the involute may be developed as follows: Let R^ 
be the radius of base circle. Let r and 0 be the radial and angular 
coordinates of a point on the involute, and let 0 be the angle displaced 
by the tangent line as shown in Figure II. 3. Hence, we have 


0 


B - ip = 6 - tan 


Jr 2 - Rh 2 

x — 


(II. 1) 


where tp is the difference between 8 and 0 as shown. 

/ 2 ~ 2 ~ 

The length of the generating line / r -R ^ is also the length of 
the ci rcumference of the base circle subtended by the angle 8. Hence, 
we have 



R, 8 or 
b 


B = 



(II. 2) 
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Figure II. 3 Involute Curve Geometry 


Then from Equation ( 1 1 . 1 ) we have 



This is the polar representation of the involute curve. The 
generating line is always normal to the involute curve and its length is 
a measure of the radius of curvature of the involute curve. For 
example, length AP is the radius of curvature of involute curve at 
point P. 

Consider the action of one involute against another involute as 
shown in Figure II. 4. The point of contact between the two involutes is 
that point where the tangents to the two curves coincide. The tangents 
to both involutes are always perpendicular to their generating lines. 
Hence, the tangents to the two involutes coincide only when the 
generating line of one is a continuation of the generating line of the 
other. Therefore, the locus of points of contact between two involutes 
is the common tangent (AB) to the two base circles as shown in Figure 
II. 4. 

When an involute is revolved at a uniform rate, the length of its 
generating line (AE) changes uniformly. Similarly, the length of the 
generating line (BD) on the mating involute changes at the same uniform 
rate. The length of the common tangent (AB) to the two base circles 
remains constant. Hence, the relative rate of motion depends only upon 
the relative sizes of two base circles. The contact takes place only 
along the common tangent. The relative rates of rotation are 
independent of the distance between the centers of the two base circles. 
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The relative rates of rotation of the two involutes may be 
represented by two plain disks which drive each other by friction. Such 
disks are known as pitch disks, while their diameters are known as pitch 


diameters. 

From Figure II. 4 we see that the intersection of the common tangent 

(AB) with the common center line (O^) of the two involutes 

establishes the pitch point P and the radii of the two pitch circles: 

R , and R„ 0 . 
pi p2 

The angle between this common center line and a line 

perpendicular to the common tangent AB is called the pressure angle. 
There is a relation between the pitch diameter, base diameter, and 
pressure angle of any given involute. That is for a given pitch 
diameter, there is a unique corresponding pressure angle: Let R bl and 
R^^ be the base radii and let R pl and R p2 be the pitch radii of first 
and second involutes. 

Let <P be the pressure angle. 

Then 


R. , = R , cos <p and R. 0 = R 0 cos 
bl pi b 2 pZ 


(II. 4) 


Next, consider Figure 1 1.5 which shows I he involute* profile change 
with roll angle. Let e represent the roll angle, and let the angular 
intervals e,, e^, and e g be equal. We note from Figure II. 5 that 
the length of the curve ab is much less than that of be; that be is 
shorter than cd; etc. Thus, when two involutes are acting against each 
other, a combined rolling and sliding action takes place due to the 
varying lengths of equal angular increments on the profiles. From 
Figure II. 5 it is evident that the rate of sliding between two involutes 
is continually varying. The sliding which is high initially reduces to 
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Figure II. 4 Action of One Involute Against Another 
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zero at the pitch point, and then changes direction, and increases 
again. 

Equations for determining the sliding speed at any point on a pair 
of involute gear teeth may be derived as follows [1]: The sliding speed 

is the difference in speeds of the ends of the generating lines of the 
involutes, projected normal to the line of action, as they pass through 
the line of action. The angular velocity of these generating lines is 
the same as the angular velocities of the gears themselves. The sliding 
speed is then determined from the products of the angular velocities and 
the lengths of the generating lines. 

Hence from Figure II. 6 the sliding speed V $ may be expressed as: 

V s ■ (R cl"l ' R C2"2> < H ' 5 > 

where R ^ and are the tooth radii of curvature of the driving and 
driven gears and where and are the angular speeds of the driving 
and driven gears. Let V be the pitchline speed. Then V may be 
expressed as 


V = R pl"l = R p2"2 


(II. 6) 


where Rp^ and Rp^ are the pitch radii of the driving and driven gears. 
Hence, by substituting for and in Equation (II. 5) the sliding 
speed becomes: 


V( 


cl 


Pi P2 


R ,R n - R ,R/-o 
=■ <(-E LSl , pl 9 


TCT, 

pl p2 


(II. 7) 
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Figure II. 5 Roll Angle and Change of Involute Profile Length 
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Figure II. 6 Involute Tooth Profile for Determining 
Sliding Velocity. 
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Next, let r^ and r^ measure the radial distance from the gear 
centers to the tooth contact point as shown in Figure II. 6. Then 
and R g 2 may be expressed in terms of r^, r^, the pressure and <p, the 
base circle radii R^ and R^ and the pitch radii as: 


R cl + R C2 = < R pl + R p2> sinf 


and 


and 


Hi 

= / r. 


R cl / r l ' R bl 


(II. 8) 


(II. 9) 


/ ? ~ ? 
R c2 ^ r 2 R b2 


= (Rpl + R p2 )sinO> - / r 1 2 3 - (II. 10) 


Finally, by substituting for R c ^ and R^ in Equation (II. 7) the 
sliding speed V g becomes (see [1]): 


V s~ R ,R 


Pi P2 


[ R p2 H 


R bl - R nl^ R pl +R p2^ n '* > 


/ 2.2 n 

■ /r r R bi ] 


ttjct <W - R pi si "^ ! 01. n) 


pi >2 


1 1 . 3 Deflection and Stiffness of Spur Gear Teeth 

The static deflection of a pair of mating teeth is assumed to be 
composed of the following components: 

1. Cantilever beam deflection of gear teeth including bending, 

shear, and compression deformation. 

2. Deflection due to rotation of tooth base because of the 
flexibility of tooth foundation. 

3. Local contact (Hertz) deformation at contact point. 
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The deflection of a gear tooth is calculated normal to the tooth 
profile. The foundation effect and the shear effect are important 
because of the "stubbiness" of a gear tooth. In this research phase we 
consider only Low Contact Ratio Gears (LCRG), with contact ratio between 

I. 0 and 2.0. However, the analysis can easily be extended to higher 
contact ratio gears. 

II. 3.1 Direction of the Applied Load on a Tooth Profile 

The application point and direction of the transmitted load by 

meshing teeth at different contact positions can be obtained through the 

i nvol utometry of the tooth profile: given the pitch radius R p , the 

pressure angle <p, and the circular pitch P c , let an X-Y coordinate 

system be defined as shown in Figure II. 7. The tooth thickness on the 

pitch circle, t , is approximately half of the circular pitch P . The 
P c 

tooth thickness t^ at any other point i on the tooth profile, with 
radius R^ from the gear center, can be expressed as 

t i = 2 R i ^p^V + inV<<> " inV0 i ^ (II. 1?) 

where 0 .j = cos'^R^/R^ ) and where inv x is the involute function defined 
as: 


i nv x = tan x - x (11.13) 

The angle x between the Y-axis and the radial line to the contact point 
is : 


t i = (tp/2R p ) + inv <J> - inv 8 i 
The angle 6 between the transmitted load at 


(II.14) 

point i and the x-axis is: 
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Figure II. 7 Direction of Transmitted Load 
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(II. lb) 


The chordal tooth thickness h at point i is 

h.j = 2R.j sim|>.j (11.16) 

Finally, the x and y coordinates at point i are: 

Y. = R i sin^ 

X. = cos (11.17) 

These coordinates will be useful in developing geometrical 
properties for calculating deflection and stiffness. 

II. 3. 2 The Tooth As a Cantilever Beam 

The involute portion of a gear rooth may be modelled as a non- 
uniform cantilever beam [9, 14]. Let l be an effective length which 
extends from the tip to the beginning of fillet area as shown in Figure 
II. 8. Let the beam be divided into a sequence of segments as shown. 
The deflection and compliance can then be obtained using the principles 
of elementary strength of materials. For each segment i, the height Y^, 
the cross-sectional area , and the area moment on inertia I ^ are taken 
as average values at both faces according to the equations: 
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Figure II. 8 Geometry and Segment Modelling of the Involute Portion 
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(11.18) 


Vi = (Vi * V, +1 ) / 2 

A, = (A i + A i+1 ) / 2 ■ 2 F V, 

I, - Hi + I H+1 ) / 2 - F(V? + Y? +1 )/3 


where F is the face width. The deflection at the loading point, in the 
direction of load, is obtained by superimposing the deflection 
contribution of the individual segments. Each segment itself is 
considered as a cantilever beam having the inner end fixed and the 
remainder of the tooth adjacent to the other end of the segment, as a 
rigid overhand. Figure II. 8 shows the applied load resolved into an 
equivalent system of forces and moments at the outer face of the 
segment. The components are: 

“id - “j cos 6 j 

“2J = V" 8 j 

M id = w j (L ij cose d • V d s1n V 


where i refers to the segment, j refers to the loading position, L. . is 

( j 

the distance from j to i, Wj is the transmitted load, Wj. and W 2 j are 
the component loads at i, and M, . is the moment at i due to the load at 

J 

j. 


II. 3. 2.1 Bending deformation 

The deformation at the load position j due to the deformation of 
segment i is the sum of the following: 
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(i) Displacements due to load kfcosB-j 


( q Bl^ij 

Wjcossj. (t 3. 

3E T. v 

e i 

(11.20) 

^ q Rl ^ i j 

^Q_s (t 2 L } 

2E P I . V i iJ 

e i 

(11.21) 


(ii) Displacements due to moments M — 


(a \ - Wjd-jfos B -j - V s i n g ■) 2 

iq B2Mj 2E I * i 

e i 


W • (L . .cosB . - Y .si nB . ) 

<Wij- - jL ^T? 1 (T ' L U> 


E e ! 1 


( 11 . 22 ) 

(11.23) 


where T.. is the thickness of segment i, q g is the displacement due to 
bending and q R is the displacement due to rotation. 

In the above expressions, E 0 is the "Effective Young's modulus of 
elasticity" whose value depends upon whether the tooth is "wide" or 
"narrow". According to Cornell [14], a wide tooth is one for which 


F/Y > 5 


(11.24) 


where Y is the tooth thickness at the pitch point. In this case the 
tooth approximates plain strain deformation and E g is 

E e = E/(l - v 2 ) (11.25) 

where E is Young's modulus of elasticity and v is Poisson's ratio. For 
a "narow" tooth, 


F/Y < 5 (11.26) 

In this case the tooth approximates plane stress deformation and E p is: 

E e - E (11.2/) 
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II. 2. 2. 2 Shear Deformation 


The shear deformation (q )., is caused by the transverse component 
of the applied load. It displaces the centerline without rotation. 
For a rectangular cross section, the shear deformation is: 

(„)..= 1.2 W jTjCOSP j^ 2.4(1+*) HjT^os Bj (I] 28) 

5 1J G A, E/ i 

II. 2. 2. 3. Axial Compression 

This axial compression (q ).. is caused by the component W, sinB. 
and is given by 

(q ) . . = si P .. Pj I i - (11.29) 

EA. 

II. 3. 2. 4 Totals 

The total displacement at the load position j, in the direction of 
the load, due to deformation of segment i can thus be found from the 
expressions : 

(1) For a wide tooth, (plane strain): 


<Vij ' < q Bl +q Rl +q B2 +q R2%tj coss j + < q ctj sin ^ 


W. [ -£2- s 

J E 


2 „ t 3 t 2, . -r . 2 


(li + T-Lii^TiLij } _ 


31, 


I . 

1 I 


cos B i sin b i T i Y i + ; 

E 


2I i 


I . 

i 


. _2 


cos ,-.?.j. ( 2 .- 4 .U + v ) T i ) + -iiQ-JLj. (li-)] (ii.3o) 


A. 

1 


*i 


24 


(2) For a narrow tooth (plane stress): 






n'ij j 


2 . 4 ( 1+ v ) + tari P -j ] (I1 . 31) 


A. 


II. 3. 3 Flexibility of the Fillet and Foundation 

The effect of the flexibility of the foundation upon the 
deformation at the load point is a function of the fillet geometry, and 
the load position and direction. This deformation is influenced by the 
effective fillet length and the fillet angle (approximately 75°). 

Figure II. 9 shows a model of a gear tooth, its fillet, and its 
foundation. Using the notation of the f i gure, def 1 ecti on in the 
direction of load at the loading point due to beam compliance of fillet 
region, q^, is given by the following expressions: 

(1) For a narrow tooth (plane stress): 


( q fb )ij w h 


( Wi 


(2) For a wide tooth (plane strain): 


2 

r COS 2 B i (T fb ) i ( Tfb>j_ t (T fh ).(L fh ) ilt (B fb )f. i _ 
J L E / T \ 




^fbh 


tanBj ^foh-Tj + Y j(L fh ) 1l + 2.4(1 + v ) + tar 2 B 


< A fb>i 


(11.32) 


( T fb>i 


(< Wij W j 


[ cos^i — ~ ^Wi^fbMWid , 2 - 4 (l^)( T f b )j 


l. E 


^ Wi 


(Tr K )?Y . 

fb'i j 




] 


cosBjSinB^ ^ 2_ 


+ (T f .) • Y • (L f .) . • sin 2 B . (T.,). n 
v fb'i y fb ; i + P J ( v _fbM ^ [ 

(i fb ) i c e ( A fbh ' J 

(11.33) 
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Figure 1 1. 9 


Model of the Fillet and Foundation 
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Similarly, the deflection due to foundation flexibility, is: 
(1) For narrow tooth (plane stress): 


L e f 


v } 


+ 1.534(1 

2.4(l + vj 


(11.34) 


(2) For a wide tooth (plane strain): 


^fe^ i j 




E F 
e 


1 -v 


tan^B- 


+ 1 . 534 ( 1 .^^] 


(11.35) 


where (L f ) i and h f are 


(L f ) i = L i + r ( s i n Yp - siny) (11.36) 

and 

h^ = h + 2r(cosy - cosy^) (11.37) 

where 

L. = (R L ) i cos(a L ) i - _ R cos a (11.38) 

A number of investigators have studied the flexibility of gear 
teeth. Experimental measurements have been made by Timoshenko and Baud 
[16], Walker [16], Buckingham [16], and Van Zandt [16]. Theoretical 
analyses have been developed by Timoshenko and Baud [16], Walker [16], 
Weber [16], Richardson [16], and Attia [16]. The results of some of 
these are compared with the results herein in Figure 11.12. It is seen 
that they have similar shapes but different magnitudes. The discrepancy 
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may be due to the following factors: 

1) Different reference points of zero deformation. 

2) Different formulas used to calculate the deformation at the 

contact area. 

3) Empirical formulas extended beyond the accuracies of 

experimental data. 

Finally, the total deflection in the direction of load due to the 
flexibility of the fillet and the foundation is obtained by adding the 
above individual deflections. That is. 


( q 2^ij ^fb^ij + ( q fe^ij 


(11.35) 


II. 3. 4 Local Compliance Due to Contact Forces 

The local contact deformation consists of two e 
Hertz (or line-contact) deformation and the other is 
of the tooth between the contact point and the tooth 
Using procedures of Lundberg and Palmgren 
compliance of cylinders in roller bearing, the local 
teeth may be approximated as: 


lements: One is the 
general compression 
centerl i ne . 

[37], for contact 
deformation of gear 


where 


, , 1.275 

'Vij ' E12 0.9 F 0.8 U 0.1 

j 


E12 = (2 El E2 ) / (El + E2) 


(11.37) 

(11.38) 


where El and E2 are the elastic modulii of the driving and driven gears 
respectively. Equation (11.37) gives the contact compliance at the jth 
calculation point. However, the compliance depends upon the magnitude 
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of the load shared at the loading position. The magnitude is unknown. 
Nevertheless, since the load in Equation (11.37) is raised to the one- 
tenth power, there is little error in replacing it by the average 
transmitted load. 

II. 3.5 Total Deflection and Stiffness 

The total deformation (q^). at load position j in the direction of 
load is the sum of all the deformations calculated above. If the number 
of segments is n, then 

<«t>j = l + + (11 - 39) 

i=l 

Figures II. 10 to 11.12 show the normalized total deformation of a 
pair of teeth along path of contact for two different cases (different 
gear radii). 

The equivalent stiffness of a gear tooth is defined as the ratio of 
the transmitted load and the total deformation: 

D 

(Kg)j = w j/ (q t )j (11.40) 

Hence the average stiffness of a meshing tooth pair is: 

< K 9>avg * 1 < K 9)j ('1.41) 


where n is the number of meshing positions. 
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Figure LI. 12 Comparison of Results for Normalized Deformation of a Tooth pair, 
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1 1 .4 Mesh Analysis 


11. 4.1 Analysis of the Meshing Cycle 

Figure 11.13 illustrates the motion of a pair of meshing teeth. 
The initial contact occurs at A, where the addendum circle of the driven 
gear intersects the line of action. As the gears rotate the point of 
contact will move along the line of action, APD. When the tooth pair 
reaches B, the recessing tooth pair disengages at D leaving only one 
tooth pair in contact. From B through P to C is then a single contact 
zone for low contact ratio gears. When the tooth pair reaches point C, 
the next tooth pair begins engagement at A and starts another cycle. 

In our analysis, the posiiton of the contact point of the gear 
teeth along the line of action is expressed in terms of roll angles of 
the driving gear tooth. This is therefore consistent with the analysis 
of gear tooth deflection and stiffness, which is also expressed in terms 
of roll angles of the driving gear tooth. 

11. 4. 2 Transmission Error, Stiffness, and Load Sharing 

Transmission error is defined as the departure of a meshing gear 

pair from the constant-angul ar speed ratio, as expected from the tooth 
number ratio. Transmission error is thus the measure of instantaneous 
variation from the ideal nominal value. Transmission error arises from 
the following sources: 

1. Combined deflection of meshing teeth 

2. Tooth spacing error 

3. Tooth profile error 

4. Runout error 
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Usually the transmission error is expressed as an angular value. 
However, linear values are convenient for use in design and analysis 
computations. Specifically, in calculating the total transmission error 
for a meshed gear pair, the linear component values for each gear are 
often the same, whereas the angular values are a spectrum, inversely 
proportional to the pitch radii. 

Transmission error and its components at the contact point may be 
expressed in the following notation: (The runout error is combined with 

the tooth spacing errors) 

E^: The deflection of gear teeth at the contact point. 

E $ : The tooth spacing error. E $ is positive if the tooth spacing 
of the driven gear is less than the base pitch P^, or if the 
tooth spacing of driving gear is greater than the base pitch 
P^. Otherwise, E s is negative. 

E p : The tooth profile error. E p is positive if material is 

removed from the surface at the contact point. Alternatively, 
E p is negative if material is added to the surface at the 
contact point. 

E^: The transmission error at the contact point. is posistive 

if the driving gear leads the driven gear. Otherwise, E^ is 
negative. 

Let the sequence of mating tooth pairs be designated by: 

a, b, c,...etc. Let the subscripts 1 and 2 designate the driving and 

driven gears respectively. Then the transmission error (E.) • and the 

t J 

load sharing W. at contact point j may be expressed as follows: For the 

double contact zone with tooth pairs a and b in contact, we have: 
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(11.42) 


= c E 3i)j + ( E d2)j- 


and 


< E t)j ■ < E Sl>J^ E 5l>j*< E 52>j + « E Sl>J +tE p2>J + < E Sl>d +(E s2>j 

/ T T 

w 


w a + w b 


(11.43) 

(11.44) 


where W is the total transmitted load. For the single contact zone, 
with only tooth pair a in contact, we have: 


< E t>j ■ 

w a = w 

J 


(11.45) 

(11.46) 


To simplify the notation, we introduce the following items: 


(E s>j 

( Vd 

< E d)j 


= UcW 


sl'j 


' $2' j 


(E P1>J + (E P2>J 
< E dl>j + < E d2>j 


Vj 


(II. 4/) 

(11.48) 

(11.49) 


where Q. is the compliance at contact point j. To simplify the 
analysis, we assume that the transmission error is the same for tooth 
pairs a and b when they are sharing the transmitted load. Hence, we 
have 


< E Vj ■ < E \>j = < E t>j <‘1. 50) 

Substituting from Equations (11.47), (11.43), (11.49) and (11.50) into 

Equations (11.42) and (11.43), we obtain the relation: 
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(11.50) 


Q a W a + (E a ) • = 0 b W b 
TJ v P ; J Vj 


(E b ) . 
v p'j 


(E b ) . 
v s'j 


Solving Equations (11.44) and (11.51) for W. and W. give the results: 

J J 


< E p)j * E < E s>j - < E p>j ] ' UI.51) 

"j ■ W> + < E p)j - < E pV - < E s>j] ' <^j b ) (11.52) 


If all errors except tooth deflection are ignored 


then W a and W b 

J J 


are: 


Wj 3 = (QjW) / ( Q j + Q b ) 

Wj b = (Q®W) / (Q a + Q b ) (11.55) 


Figure 11.14 shows typical stiffness and load sharing 
characteristics of a low contact ratio gear with tooth deflection but no 
other transmission errors. Let a series of mating tooth pairs be 
denoted as a, b, c, d and let points A, B, P, C, D be the same as those 
in Figure 11.13. Then, AB and CD represent the double contact regions, 
BC represents the single contact region, and P is the pitch point. 

The stiffness values at double contact regions are clearly much 
higher than those at single contact regions. When gears rotate at 
appreciable speed, this time-varying stiffness as shown in Figure 11.16 
is the major exictation source for the dynamic response of gear systems. 

The total transmitted load is shared between two pairs of teeth 
within double contact regions. The magnitude of the load shared by each 
individual tooth pair depends on the stiffness of that tooth pair. The 
higher the stiffness, the higher the shared load. 
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II. 5 Lubrication and Friction 


II. 5.1 Boundary, Mixed, and Hydrodynamic Lubrication 

A purpose of lubrication is the reduction of friction between mating 
surfaces. The relative motion between mating gear teeth produces 
sliding and rolling friction which can be reduced but not totally 
eliminated. The friction developed depends upon the following: 

(1) The characteristics of the lubricant. 

(2) The pressure exerted between the contacting bodies. 

(3) The nature and condition of the contacting surfaces. 

(4) The ambient and surface temperatures . 

There are four basic cases of lubrication which may be described as 
fol lows: 

1. Boundary Lubrication : In this case the friction is mainly due to 

the interaction of asperities of the contacting surfaces. 

2. Mixed Lubrication [Partial Elastohydrodynamic(EHD) Lubrication]: In 

this case the friction is created by both the interaction of 
asperities and by elastohydrodynamic effects. 

3. Elastohydrodynamic (EHD) Lubrication : In this case there is no 

interaction of the asperities. The friction is a function of the 
fluid properties of the lubricant and in certain situations, the 
surface quality. The mating surfaces deform elastically. 

4. Hydrodynamic Lubrication : This case is the same as EHD Lubrication 

except that there is no deformation of the mating surfaces. 

Figure 11.15 illustrates that these four cases, and Figure 11.16 shows a 
representation of the friction coefficient for the four cases. 
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Figure 11.15 Lubrication Cases: (1) Hydrodynamic ; 

(2) Elastohydrodynamic (EHD); (3) Mixed; 

(4) Boundary. 



11.16 Coefficient of Friction of Different Lubrication Cases 
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An important parameter for separating regimes is tjpe specific film 


thickness X, defined as: 


X 


h . / S" 

min 


(11.56) 


where h m ^ n is the EHD minimum oil film thickness and S" is the composite 
surface roughness. S" may be expressed as: 

S" = A\ + S Z 2 (11.57) 

The surface roughness varies for S 1 and are normally 
arithmetic-average (AA) values. [Formerly, root-mean-square (RMS) 
values were used for surface finish. The RMS value is approxi matl ey 
1.11 times the AA value.] 

For Boundary Lubrication , xis less than 1 (x < 1 ) , and the gears may 
be thought of as running wet with oil. In this case the thickness of 
the oil film is quite small compared with the surface roughness. 
Essentially full surface-to-surface contact is obtained in the Hertzian 
contact band area. This is typical of slow-speed, high-load gears 
running with a rough surface finish. Hand-operated gears in winches, 
food presses, and jacking devices are typical of such gears. 

For Mixed Lubrication x is between 1 and 5 (1<X<5), and partial 
surface-to-surface contact. The asperities of the tooth surfaces 
collide with each other, but substantial areas of the surfaces are 
separated by a thin film. Mixed Lubrication is typical of medium-speed 
gears with a good surface finish, which are heavily loaded, and which 
are running with relatively viscous oil. Most vehicle gears run under 
mixed lubrication. 
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For EHD Lubrication X is between 1 and 10 (1 <A<10) and for 


Hydrodynamic Lubrication X is greater than 10 ( X > 10). In these cases, 
the oil film between contacting tooth surfaces is thick enough to avoid 
surface-to-surface contact. Even high asperities miss each other. 
Precision, high-speed gears generally operate under these lubrication 
conditions. Turbine-gears in ship drives, gears in electric generators, 
and gears in compressors are examples. In the aerospace industry, 
turboprop drives and helicopter main rotor gears operate in this region. 
An exception may be some final-stage gears which may run slowly enough 
to be in the mixed lubrication region. 

References [17, 18, 19 and 31] provide additional details about the 
lubrication cases and about surface finishes. 

II. 5. 2 Coefficient of Friction, Frictional Torques 

There is a lack of agreement as to the proper form and variation of 
the friction coefficient. Buckingham [1] has recorded a semiempirical 
formula for the friction coefficient f as: 

f = 0.05e' 0,125V s + 0.002 /T (11.58) 

where is the sliding speed measured in in/sec. An analogous 
expressions has been developed by Benedict and Kelly [22] and by 
Anderson and Lowenthal [23]: 

f = 0.0127 log (45.94 W/F Mq V Vj[) (11.59) 

where W is the applied load, measured in lb, F is the face width, 
measured in inches, i$ the rolling velocity, measured in i n/ sec , and 
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Figure 11.19 Friction Torques Variation Along the Contact Path 
Using Buckingham's Model 


45 













is the lubricant viscoscity, measured in lb sec/in. 

Figures 11.17 and 11.18 show graphs of the friction coefficient as 
given by Equations (11.58) and (11.59) as a function of the roll angle. 
Figures 11.19 and 11.20 show the resulting effect upon the friction 
torque. 

Finally, it should be noted that Equation (11.59) is valid only 
when the specific film thickness A [See Equation (11.56)] has values 
between 1 and 2. However, this is the range of many gear applications 
of interest. Moreover, extending the use of Equation (11.59) to systems 
with higher values of A should not introduce serious errors. 

II .6 Shafts and Connected Masses 

The shafts and connected masses are also components of the gear 

system. The parameters related to dynamic system response are the polar 

mass moment of inertia, and the stiffness of shafts and masses. 

For bodies with annular cross-section, outside diameter D , inside 

o 

diameter , the polar area moment of inertia J is: 

J = 52 (°o ' (II - 60) 

The torsional shaft stiffness K is then given by: 

K = JG/Jt (11.61) 

where G is the shear modulus and l is the shaft length. 

The polar mass moment of inertia I is given by the expression: 

I = L m fn 2 - n 2 l (11.62) 

a ■' o i' ' ' 

where m is the shaft mass. 
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III. SYSTEM DYNAMIC ANALYSIS 


111.1 Dynamic Analysis and Computational Procedure 

111. 1.1 Mathematical Modelling and Equations of Motion 

The system model as shown in Figure II. 2 is used for the dynamic 
analysis. Its mathematical model, shown in Figure II. 1, and the 
resulting governing equations are developed using basic gear geometry 
and elementary vibration principles. The equations may be expressed as 
follows: 

J M 0 M + C s1^M " ®1^ + K sl^ 9 M ' 6 1 ) = T M 

J l' 0 l + C sl^l ' 9 M^ + K sl (e l ” e M^ + C g^^ R bl 9 l ‘ R b2®2^ 

+ ^ R bl^ R bl 9 l ' R b2 6 2^ = Tfi^) 

J 2 § 2 + C s2^2 - V + K s2< 6 2 ' V + C g (t)[R b2®2 ' R bl 5 ^ 

+ Kg^ *- R b2 ^ R b2 6 2 " R bl 9 l^ = T f2^ 


J L § L + C s2 (0 L " 6 2 } 


+ K s 2 (e L - e 2 ) - - t l 


where the notation is 


(III.l) 


9, 0, 0 angular displacement, velocity and accelera- 

tion 

e M’ 0 i* 0 2’ 0 L angu 1 a r displacement of motor, gears 1 and 2, 

and the load 
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J M’ J l’ J 2’ J L mass moments °f inertia of the motor, gears 

1 and 2 and the load. 

C s r C S 2 > C ( t ) damping coefficients of the shafts and gear 

teeth . 

T M » T ^ , T ^ j ( t ) , T ^ ( t ) input and output torques and 

frictional torques on the driving (1) 
and driving (2) gears. 

The procedures developed in Part II are used to determine the 
spring stiffness, modelling the elastic behavior of the gear teeth. The 
equations of motion will thus be dependent upon excitation terms due to 
transmission errors. However, it is a difficult task even under ideal 
geometry conditions to obtain a precise modelling of the dynamic 
loading. The vibration of the gear tooth is affected by irregularity in 
motion due to manufacturing errors and by deflections of the tooth due 
to the applied loads. The motion is further affected by other teeth 
leaving contact. The motion becomes more complicated when the teeth 
begin to assume the entire load, resulting in a complex vibratory motion 
taking place in the zone of single contact. 

The meshing process leads to instantaneous load fluctuations in the 
teeth even under constant loading conditions. Also, the magnitude of 
the load and the fluctuations are influenced by the damping effect of 
the lubricant, and the proximity to system natural frequencies of the 
operating frequencies. 

Several assumptions are needed to simplify the analysis: 

1. The dynamic process is defined in the rotating plane of the 
gears. The out-of-plane twisting and misalignment are 
neglected. 
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2. Damping due to lubrication of the gears and shafts are 

expressed as constant damping coefficients. 

3. The differential equations of motion are developed using the 

theoretical line of action. 

4. Gear tooth stiffness is evaluated according to dynamic 

considerations: That is, the reference point for tooth 

deflection is assumed to be located at the radius of gyration. 

5. Low contact ratio gear pairs are used in the analysis. The 

contact ratio is taken between 1 to 2. 

III. 1.2 Computational Procedure 

Equations ( 1 1 1 . 1 ) are solved using a 1 inearized-iterative proce- 
dure. The linearized equations are obtai ned by dividing the mesh period 
into n equal intervals. In the analysis a constant input torque T M is 
assumed. However, the output torque T L is assumed to be fluctuating due 
to damping in the gear mesh, friction, and time-varying mesh stiffness. 
Finally, a true involute profile is assumed. 

Initial values of the angular displacements are obtained by 
preloading the input shaft with the nominal torque carried by the 
system. Initial values of the angular speeds are taken from the nominal 
operating speed of the system. 

The iterative process is as follows: The calculated values of the 

angular displacements and angular speeds after one period, are compared 
with the assumed initial values. Unless the differences between them 
are sufficiently small, the procedure is repeated using the average of 
the initial and calculated values as new initial values. The criteria 
for repeating is: 
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and 


|e_ - el > 0.001 0 

'no' o 


1 6 - ej >0.001 e 

'no' o 


(III. 2) 


The term (R bl e 1 - R^ ® 2 ) ' in the equations of motion represents 
the relative dynamic displacement of the gears. Let Bh represent the 
backlash. Let gear 1 be the dri vi ng gear. The f ol 1 owi ng condi tions 
can occur: 

Case (i) 


R bl 9 1 ' R b2 0 2 > 0 


(III. 3) 


This is the normal operating case and the dynamic mesh force is 
defined as: 


W dl = V^^bl 6 1 ‘ R b2 0 2^ + C g^ t ^ R bl 0 l ‘ R b2®2^ 
and 

W d2 = W dl 
Case ( ii ) 

R bi e i " R b2 0 2 ^ ^ an< ^ I R b 1 ® 1 ~ R b2 0 2 I ^ 

In this case, the gears will separate and the contact between the 
gears will be lost. Hence, 

W d = 0 (III. 7) 

Case ( i i i ) 

R bl 0 l ' R b2 9 2 < 0 and l R bi 0 i ' R b 2 6 2 1 > Bh (HI. 8) 


(111. 4) 

(111. 5) 

(111. 6) 
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In this case, gear 2 will collide with gear 1 on the backside. Then, 


W dl = K g ( t ) c ( R b2 0 2 " R bl 6 1 ^ " Bt1 ^ + C g^ t ^ R bl®l ‘ R b2*M 

(HI. 9) 

and 

“d2* H dl (III ' 10 > 

The equations of motion contain damping terms for all components in 
the system. The damping in the shafts is due to material damping. 
Experiments have suggested that this damping is between 0.5% and 0.75 of 
the critical damping (See Reference [35].) Hence, the effective damping 
of shafts is taken as 

c si = 2 £ s / K s1 /(1/J d + 1/J 1 ) (III. 11) 

and 

C s2 = 2C S /k s2 /(l/J L + 1/J 2 ) (III. 12) 


where represents the critical damping ratio of shafts with value: 
0.005. 

Similarly, the effective damping of the gear mesh is taken as: 


C g = 26[ Y blRb f JlJ ^ 

* R bl J l +R b2 J 2 


(III. 13) 


where S represents the critical damping ratio. Measurements of geared 

systems show £ to range between 0.03 and 0.17 (See References [10, 11]). 

In Equation (III. 13) the average gear mesh stiffness is used and the 

equivalent masses of the gears are concentrated at the base circle to 

reflect their effects along the line of action. 
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The flow chart of the computational procedure is shown in Figure 
III. 2. 

III. 2 Natural Frequencies 

The natural frequencies of the system may be obtained by examining 
the undamped equation of motion written in matrix form: 


[J][e] + [K][e] = [0] 


(III. 14) 


where the inertia matrix [J] is 


[J] - 


0 

0 

0 


0 

J ] 

0 

0 


0 

0 

J, 


0 

0 

0 

J, 


(III. 15) 


and the stiffness matrix [K] is 


[K] = 


v sl 




■ K sl K sl +(K g ) avg R bl '(V avg R bl R b2 0 


0 -( K g)avg R bl R b2 K s2 + ( K g)avg R b2 ' K s2 


's2 


v s2 
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The average value of gear mesh stiffness is taken as the sum of the 
discrete tooth stiffness values over one mesh cycle divided by the 
number of mesh positions in the cycle. 

Equations (III. 14) may thus be used to determine the eigenvalues 
and thus the natural frequencies of the system. 

III. 3 Parametric Study 

A brief parametric study was conducted to determine the effects of 
the rotating speed, the load, the stiffnesses, the inertias, and the 
contact ratio on the dynamic load. The system examined consisted of two 
identical gears having the following properties: 

Number of teeth: 36 

Diametral pitch: 8 

Pitch diameter: 4.5 in 

Pressure angle: 20° 

Face Width: 1.0 in 

Moment of inertia: 0.02947 in lb sec^ 

Stiffness: 3.5355 x 10 6 lb/in rad 

Damping ratio: 0.10 

The shaft stiffness, inertias were: 

Shaft stiffness: 10081 in Ib/rad 

o 

Motor and load inertia: 0.08841 in lb sec (each) 

Finally, the first three system critical speeds were: 
o) : 89.6 rpm, 179.1 rpm, and 8688 rpm 
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For gears with true involute profile under normal operating 
conditions, the principal source of dynamic excitation is the time 
varying tooth stiffness due to alternating single and double contact of 
the tooth pairs. The gear system vibration depends upon the frequency 
of this parametric excitation, and thus also the operating speed. 

III. 3.1 Dynamic Load as a Function of Roll Angle 

Figure III. 3 shows the variation of dynamic load response for a 
pair of teeth as a function of roll angle. At lower speeds where the 
excitation frequency is much lower than the resonance frequency, the 
dynamic load response is basically a static load sharing in phase with 
the stiffness change, superimposed with an oscillatory load at a 
frequency corresponding to the natural frequency. 

At higher speeds, close to the natural frequency, the dynamic load 
variation becomes so abrupt that it produces tooth separation. The peak 
dynamic load is much higher than the static load and is very likely a 
source of gear noise and early surface fatigue. 

Figure III. 4 shows the dynamic load response as a function of 
operating speed. The major peak occurs at the natural frequency. The 
second major peak occurs at the principal parametric resonance frequency 
due to the parametric excitation of the time varying tooth stiffness. 
This parametric resonance frequency is about half the system natural 
frequency [27]. Finally, the dynamic factor is defined as the ratio of 
dynamic load to static load. 

Observe that for speeds above the natural frequency, the dynamic 
response decreases steadily in the same manner as with elementary 
vibratory systems. 
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Figure III. 5 shows a three dimensional representation of the system 
dynamic response. The horizontal axis represents the operating speed, 
and the contact position along the tooth profile. The total number of 
contact positions is 121. The vertical axis is the dynamic factor, 
defined by the ratio of the dynamic load to the static load. 

Figure III. 6 presents a contour plot of the system dynamic 
response. 

The shaded areas, represent regions where tooth separation occurs. 
They are located in the double contact regions. At near resonance 
speeds the vibration motion of gear bodies exceeds the deflection of the 
meshing tooth pairs, inducing tooth separation. 

As the speed increases, the dynamic response also shows a phase 
shift towards the higher numbered contact positions. This phenonemon 
can be seen by noting that the highest dynamic effect (at speeds ranging 
from 600 rpm to 10400 rpm) occurs at different contact positions 
beginning from position 51 and gradually changing to position 85. This 
is an important factor when computing the root stresses. 

111. 3. 2 Effect of the Applied Load 

Figure III. 7 shows the maximum dynamic load on the system as a 
function of the transmitted loads for different speeds. 

The deflection of gear teeth changes the tooth stiffness and load 
sharing between the teeth. The separation of gear teeth due to 
vibration is resisted by the transmitted load. Therefore, the dynamic 
load effect is decreased as the applied load is increased. 

111. 3. 3 Effect of Damping 

The damping due to viscous friction governs the dynamic load 

response of the gear system. Figure III. 8 shows the effect of damping 

on the dynamic load response. It is seen that the damping has a greater 
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Figure III. 7 Effect of Applied Load on the Dynamic Load Response 
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Figure III. 8 Effect of Damping on the Dynamic Load Response 
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Figure III. 9 Effect of Tooth Stiffness on Dynamic Load Response 
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Figure III. 10 Effect of Gear Inertia on the Dynamic Load Response 


66 



effect near the reasonance frequencies than at frequencies away from re- 
sonance. 


111. 3. 4 Effect of Stiffness and Mass Moment of Inertia 

A. Stiffness : Changes in shaft stiffness have a minor effect on the 
system dynamic response. However changes of tooth stiffness have a 
major effect on the response. Figure III. 9 shows that the higher the 
tooth stiffness the lower the dynamic response. This is consistent with 
observations that as the tooth stiffness increases the mass effect of 
gear bodies on the system dynamics is reduced. Figure III. 9 also shows 
that system resonance frequencies are increased as the tooth stiffness 
increases. This is a potentially useful effect for the design of gear 
systems. 

B. Mass Moment of Inertia : The effect of the shaft masses is small 
compared to that of the gears. Figure III. 10 shows that as the gear 
inertia is decreased the dynamic response is also decreased. 

111. 3. 5 Effect of Diametral Pitch and Contact Ratio 

For gears with different diametral pitches, the dynamic response is 
different due to the change in contact ratio. Gears with a finer pitch 
have a higher contact ratio. Since the contact ratio is a measure of 
the duration of the load being shared by more than one pair of teeth, it 
has a significant effect on the system dynamic response. 


Figure III. 11 shows a comparison between gears having different 
diametral pitches. The finer pitch gears, having a higher contact 
ratio, have a smaller dynamic load effect than the coarser pitch gears. 
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Figure III. 11 Effect of Diametral Pitch on Dynamic Load Response 
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III. 4 Comparison with Experimental and Empirical Data 

At an American Gear Manufacturer's Association (AGMA) meeting in 
1927 A. A. Ross [24] introduced the following empirical formula for the 
dynamic factor K y : 


78 

78 + /T 


(III. 17) 


where v is the pitch line speed measured in ft/min. This expression 
came to receive acceptance as a standard factor used by the AGMA. 

A similar factor for use with higher precision gears was introduced 
by Wellauer [25] in 1959: 


K v 2 — 

v 78 + iT 


(III. 18) 


Equation (III. 17) and (III. 18) are recognized as being conservative 
when applied with very high precision gears. That is they predict 
dynamic laods which are larger than the physical loads. 

Buckingham [1] also developed an expression for the dynamic load in 
terms of the pitch line speed, the applied load and other geometrical 
and physical factors. Buckingham's formula is 


W 


d 


■ " + t f a< 2f 2 



(III. 19) 


where W . is the dynamic load, W is the applied load and the factors f 

Q a 

and are 
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and 


f a = W< f ! + f 2) 


(III. 20) 


f 2 = FC + W (III. 21) 

where F is the face width in inches and C is a deformation factor given 
by: 


C = 0.111 e E 1 E 2 /(E 1 + E 2 ) (111.22) 

where e is a profile error factor and and are elastic contacts of 
the gears. The factor f, is 

fj = 0.00025[(R 1 + R 2 )/R 1 R 2 ]mV 2 (III. 23) 

where Rj and R 2 are pitch radii of the gears, measured in inches, m is 
the effective mass of the gears and, as before, V is the pitch line 
speed in ft/mi n. 

Finally, Kubo [26] has recorded results of extensive experimental 
studies on dynamic effects for high-precision spur gear systems. 

Figure III. 12 shows a comparison of the AGMA high-precision formula 
[Equation (III. 18)], Buckingham's formula, Kubo's results, and the 
results of the computer simulation. 
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Figure III. 12 Comparison of Experimental, Empirical, and 
Computer Simulation Results 
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IV. DISCUSSION 


In reflecting upon the analysis there are several remarks which 
might summarize the findings: 

1. The model of Figure II. 2, Chapter 2 is relatively simple. 
Indeed, it consists of only two spur gears and two shafts. Still, many 
important aspects of dynamic loading are manifested by the model. 

2. The model includes the effects of: 

a. Masses and movements of inertia of the gears and shafts. 

b. The stiffnesses of the gears and shafts. 

c. The operating speeds. 

d. The friction and damping. 

e. The tooth geometry. 

3. The tooth stiffness is based upon a finite segment modelling. 
It is developed by using fundamental principles of mechanics of 
materials. It is dependent upon the tooth geometry--that is, the 
involute form. 

4. The loading is found to be significantly affected by the 
contact ratio. 

5. Tooth separati on--l eadi ng to impact--occurs in the double 
contact region since the deflections are smallest in that region. 

6. Dynamic load generally increases with operating speed until a 
system natural frequency is reached. The dynamic load decreases rapidly 
beyond the natural frequency. 

7. The dynamic factor is largest for contact points near the tooth 

tip. 
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8. The applied load has a minor effect upon the dynamic factor. 
However, the dynamic factor decreases with applied load. 

9. Damping and friction decrease the dynamic factor with the most 
dramatic effects occuring near the system natural frequencies. 

10. Tooth stiffness has a significant effect upon the dynamic 

factor: The higher the stiffness, the lower the dynamic factor. Also, 

the greater the stiffness, the higher the rotating speed of peak 
response. 

11. Shaft moment of inertia has a minimal effect upon the dynamic 
factor. However, the greater the inertia, the greater the dynamic 
effect. 

12. Finally, for increased dimetral pitch, that is, high contact 
ratio gears, the dynamic factor will be decreased. 
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